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Abstract. Parameters of magnetic flux distribution inside low-latitude coronal 
holes (CHs) were analyzed. A statistical study of 44 CHs based on Solar and 
Heliospheric Observatory (SOH[0)/MDI full disk magnetograms and SOHO/EIT 
284A images showed that the density of the net magnetic flux, Bnct, does not 
correlate with the associated solar wind speeds, Vx- Both the area and net flux 
of CHs correlate with the solar wind speed and the corresponding spatial Pearson 
correlation coefflcients are 0.75 and 0.71, respectively. A possible explanation for 
the low correlation between Bnet and is proposed. The observed non-correlation 
might be rooted in the structural complexity of the magnetic field. As a measure 
of complexity of the magnetic field, the filling factor, /(r), was calculated as a 
function of spatial scales. In CHs, /(r) was found to be nearly constant at scales 
above 2 Mm, which indicates a monofractal structural organization and smooth 
temporal evolution. The magnitude of the filling factor is 0.04 from the Hinode 
SOT/SP data and 0.07 from the MDI/HR data. The Hinode data show that at scales 
smaller than 2 Mm, the filling factor decreases rapidly, which means a mutlifractal 
structure and highly intermittent, burst-like energy release regime. The absence 
of necessary complexity in CH magnetic fields at scales above 2 Mm seems to be 
the most plausible reason why the net magnetic flux density does not seem to be 
related to the solar wind speed: the energy release dynamics, needed for solar wind 
acceleration, appears to occur at small scales below 1 Mm. 
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1. Introduction 

Since tiie Skylab mission in early 1970s, it is believed that the coronal 
holes seen on the surface of the Sun are related to the high-speed 
streams of the solar wind and they might be a cause of geomagnetic 
disturbances (see, e.g., Sheeley et al. (1976) and references therein). A 
possibility of ground-based observations of coronal holes in the spectral 
line He 1083.0 nm (Harvey et al., 1975) stimulated the interest to 
the problem. Numerous sophisticated models were proposed to explain 
the CHs formation and evolution {e.g., Wang and Sheeley, 1991; Fisk, 
1996, 2001, 2005; Fisk et al, 1999; Wang et al, 2000; Schrijver, 2001; 
Schrijver and Title, 2001; Schrijver et al, 2002; Schrijver and DeRosa, 
2003). 

© 2009 Kluwer Academic Publishers. Printed in the Netherlands. 



msl4.tex; 18/08/2009; 18:01; p.l 



2 



During the last decade, a connection between the coronal holes and 
various heliospheric phenomena, such as, high speed streams, corotat- 
ing interaction regions, long-living geomagnetic storms without CMEs 
(see, e.g., Vrsnak et ai, 2007 for references), became well established 
and stimulated elaboration of approaches to forecast the geomagnetic 
response to the transit of a coronal hole over the solar disk. An area 
occupied by coronal holes turned to be a very fruitful parameter. Rob- 
bins et al. (2006) suggested an empirical model to predict the solar 
wind speed at 1 AU from measurements of the fractional area occupied 
by a CH inside a 14° - sectoral region centered at the central meridian. 
A similar technique was independently applied later by Vrsnak et al. 
(2007) and further extended to forecast, along with the solar wind 
speed, other parameters of the solar wind, such as proton density, 
temperature and magnetic field strength at 1 AU. 

Several groups analyzed properties of the magnetic fields inside 
coronal holes {e.g., Harvey and Sheeley, 1979; Harvey et al, 1982; 
Obridko and Shelting, 1989; Bumba et al, 1995; Obridko et al, 2000; 
Wiegelmann and Solanki, 2004; Wiegelmann et al, 2005; Abramenko 
et al, 2006; Hagenaar et al, 2008). The most comprehensive study 
made so far on the net magnetic fluxes and averaged flux densities in 
coronal holes was presented by Harvey and colleagues (Harvey et al., 
1982), which was based on the Kitt Peak full disk magnetograms (1 
arcsec pixel size) and hand-drawn maps of coronal holes derived from 
the He 1083.0 nm data. These authors analyzed the ascending phase of 
the 21st cycle and reported an increase in the averaged net flux density 
as the solar activity intensified. This was explained by an extra flux 
deposited into low-latitude coronal holes by decaying active regions. 
Since then, the density of the net flux inside CHs was adopted as a 
representative characteristic of the magnetic flled inside coronal holes. 

As long as it is the magnetic field that is ultimately responsible for 
energetics in a coronal hole, it would be interesting to explore how 
the density of the net magnetic flux is related to the speed of the 
fast solar wind. We performed such a statistical study on the basis 
of the MDI full disk magnetograms (Section 2) and found a rather 
surprising result: the density of the net flux is not correlated with 
the solar wind speed measured at 1 AU. What could be a reason for 
that? We suggested that the reason might be that the density of the 
magnetic flux, derived from low-resolution data and averaged over the 
CH's area, is not a suitable parameter to qualify energetics in CHs. The 
ultimate reason for the observed non-correlation might stem from the 
structural organization of the magnetic flux and its multifractal nature 
at small scales. This encouraged us to study multifractal properties of 
the magnetic flux inside a coronal hole taking advantage of the high 
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resolution Hinode observations of the magnetic field (Section 3). Our 
final section represents summary and discussion of the results. 



2. Averaged Net Flux Density Versus the Solar Wind 

Speed 

2.1. Event Selection 

Here we focus on the distribution of the magnetic flux density of 44 
CHs observed between March 2001 and July 2006 at low solar latitudes 
near the center of the solar disk (Table I). To avoid influence of the 
projection effect, we required that the angular distance, 9, from the disk 
center to the center of gravity of a CH should not exceed 20 degrees. 
The values of 9 (positive when the gravity center was located in the 
northern hemisphere) are shown in the 8-th column of Table I. In the 
process of event selection we discarded all CHs that had more than 5% 
of their area outside a circle of 30 degree radius centered at the solar 
disk center. 

We also made sure that the solar wind speed profiles were not 

contaminated by a possible influence of ICMEs. Following Arge et al. 
(2004), we selected only the events when, at the arrival time of the 
streams, the value of the plasma beta 

{http : / /omniweb.g.sf c.nasa. gov/ form/dxl. html) was well above 0.1. 

For the selected 44 CHs (Table I), we utilized the following data 
sets: i) Michelson Doppler Imager (SOHO/MDI; Scherrer et ai, 1995) 
1 min averaged full disk magnetograms (spatial resolution 4 arcsec and 
the pixel size of 2 arcsec); ii) Fc XV 284 A images from the EUV 
Imaging Telescope (SOHO/EIT; Delaboudiniere et al, 1995), and iii) 
ACE/SWEPAM and MAG measurements of the solar wind at 1 AU. 

Each CH was determined in an EIT Fe XV 284 A image as an 
area with pixel intensity below a certain threshold, namely, 80 DN. 
The threshold level was first determined by trial-and-error method 
for one CH and it remained the same for all the CHs analyzed here. 
Next, we co-aligned the EIT image, taken at the time when the CH 
passed the central meridian with the closest MDI magnetogram. The 
CH boundary was then mapped on that magnetogram, and calculations 
of the magnetic field statistics were performed for an ensemble of pixels 
enclosed by this boundary (Figure 1). In cases when the time difference 
between the EIT image and the corresponding magnetogram exceeded 
1/2 hour, we aligned the images taking into account the differential 
rotation of the Sun. We discarded all events where the time difference 
between the EIT and MDI images exceeded 4 hours. 
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One more criteria was used to ensure that a given CH is indeed 
associated with the feature observed in the solar wind speed profile. In 
case of the positive association, the polarity of the Bx component of 
the solar wind magnetic field, measured in the GSE coordinate system 
by ACE/MAG, should be opposite to that of the open flux in the base 
of the CH because the x-axis in the GSE system points from the Earth 
toward the Sun, while the positive magnetic field on the solar surface 
is directed outward from the Sun. Only those CHs, for which magnetic 
polarities satisfied the above condition, were included into this study. 




Figure 1. SOHO/MDI magnctogram (left) and the corresponding SOHO/EIT 284 
Aimage taken at 19:06 UT on 2 March 2002 (right). The contour indicates the 
boundary of coronal hole CH201 (see Table I). 



2.2. Solar Wind Speed Derivation 

For each CH in our data set, we determined the solar wind speed, Vx- 
We utilized the 64 s averaged time profiles measured in the GSE coor- 
dinate system with ACE/SWEPAM instruments. Note that the solar 
wind speed acquires negative values in the GSE coordinate system. 

The arrival time at 1 AU of the solar wind associated with the CH 
was determined as the moment t, when a function 

Ar(t) = Vx{t - to) - lAU (1) 

turns into zero. Here, (3rd column of Table I) is the CH culmination 
time, and t = tx (4-th column of Table I) is the arrival time. 

We then accepted that the association between the CH and the 
solar wind feature is reliable when calculated from Equation (1), 
falls into the well pronounced "dip" in the observed solar wind speed 
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profile. Figure 2 shows an observed time profile of the solar wind speed 
and the corresponding profile of Ar(t) associated with the CH shown 
in Figure 1. 




Day of year 2002 

Figure 2. Observed time profile of tiie Vx component of the solar wind speed (red 
curve) as measured by ACE satellite from OOOOUT on 28 February 2002 until 2359UT 
on 7 March 2002. The black curve represents a smoothed over 21-point observed 
speed profile. The vertical dotted line indicates the CH culmination time, to, while 
the vertical solid line indicates the solar wind arrival time, Ia- The arrival time was 
determined as the moment when parameter Ar (double green curve) is equal to zero. 
The blue thick horizontal line segment shows the time interval used to determine the 
average, for a given event, solar wind speed, (14) (indicated by the vertical position 
of the blue line segment). 

The magnitude of the solar wind speed, (Kc), was determined by 
averaging the observed time profile of Vx{t) over a time interval that 
includes the minimum of the observed speed profile. This time interval 
(marked by the thick blue horizontal segment in Figure 2) was defined 
in the following way. First, we applied a 21-point running averaging to 
the observed Vx{t) profile (red curve in Figure 2). This smoothed profile 
(black curve) was then used to determine i) the level of undisturbed 
slow-wind speed, Vq, preceding the CH, and ii) the maximum speed 
inside the CH, l^nax- To obtain the averaged solar wind speed, (14), 
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we averaged all data points inside an interval where Vx (t) is lower than 
Vmax+ (l^axl ~ |M3|)/4. This threshold was chosen by a trial- and-crror 
method. The resulting averaged solar wind speed is indicated in Figure 
2 by the vertical position of the blue line segment. We would like to 
emphasize here that we did not include in our study those CHs whose 
associated speed profile was so complex that no well defined minimum 
could be chosen. 

For each CH in our data set, the area of the CH (in pixels of the 
MDI full disk magnetogram) and the associated solar wind speed are 
listed in the 7-th and 9-th columns of Tabic I, respectively. 

The above routine allowed us to determine the arrival time and the 
magnitude of the solar wind speed individually for each coronal hole. 
This technique differs from that applied by Robbins et al. (2006) and 
Vrsnak et al. (2007), where the mean arrival time, determined from the 
cross-correlation technique, was utilized. 

2.3. Magnetic Parameters 

The above selection routine left us with a set of 44 CHs. MDI full 
disk 1 minute average magnetograms were available for all of them. 
However, the full disk data of 1 minute cadence, as well as the obser- 
vations in the MDI high resolution mode, were not available for all of 
them. Since we are interested in measurements of the net magnetic flux 
density, we decided first to focus on whether the measurements of this 
parameter from 1 minute full disk magnetograms are reliable. To do 
that, we analyzed coronal hole CH 204 from our data list (see Table 
I) , for which nearly simultaneous observations in two MDI modes were 
available: 1 minute cadence full disk and high resolution modes. Inside 
this large coronal hole, located at the solar disk center, we outlined a 
rectangular area of n x m pixels and calculated magnetic parameters 
for five different magnetograms. Results are compiled in Table II, where 
the first row is the result from one original individual MDI full disk 1 
minute magnetogram. The second row represents results from the same 
magnetogram, but smoothed in both directions with a boxcar average 
of three pixels. The third row are parameters from a magnetogram 
derived from averaging of five consecutive original 1 min magnetograms, 
whereas its smoothed version is presented in the next row. The last row 
represents results from the high resolution MDI magnetogram recorded 
at the same place on the Sun with pixel size of 0.6 arcsec. 

From each of the above magnetograms, we first calculated the net 
flux density, Bnet = I] -B||/(n x m) (2nd column in Table II), inside the 
box. We note that the net flux is obtained by summing all pixel values 
with their sign. This procedure results in cancellation of the bulk of 
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sign- variable noise. Magnetic flux confined in loops, closed within the 
CH, will also be canceled. This allowed us to assume that in CHs, the 
net flux density also represents the open magnetic flux density. 

The density of the total unsigned flux, iJ^ot = ^ "^)) 

is shown in the 3rd column of Table II. The 4-th column represents 
the density, B^, of the noise in magnetograms derived from the power 
spectrum calculations as it was suggested by Longcope and Parnell 
(2008). To derive this parameter, we calculated a ID power spectrum, 
S (k) , from each magnetogram and then calculated the noise density as 
S„ = (7rA;25'(/cJ)i/2^ where k,, = 271^/"^ Ax is the maximum wave num- 
ber. The imbalance of the magnetic flux is shown in the 5-th column. 
The 6-th column represents the density of the net flux calculated over 
the area where |-B||| < SB^ (we denote SB^ = p), and the last column 
shows the fraction of the box area, where the flux density exceeds the 
triple noise level, > p. 

Data of Table II indicate that the density of the net flux, B^et , does 
not vary much (by less than 2%) with temporal/spatial averaging and 
does not depend on the noise level. At the same time, the density of 
the unsigned flux, B^ot, depends signiflcantly on the noise level. (Note 
that the magnitude of B^ is in good agreement with the results from 
MDI full disk data noise reduction reported by Hagenaar et al. (2008)). 
Last two columns in Table II show that the bulk of the magnetic flux 
in CHs is concentrated within a small (a few percents) fraction of the 
CH's area. Vast zones of low fluctuations contribute only about 1 G 
into the resulting magnitude (4 G) of the net flux density. 

This experiment shows, first, that an averaging procedure is undesir- 
able when one intends to analyze fine structures of the magnetic field. 
Second, that the net fiux density in coronal holes is measured with the 
same level of confidence from the original and averaged magnetograms. 
High noise level in original magnetograms does not influence much on 
the i?nct calculations. The reasons for that arc: i) high imbalance of the 
magnetic flux inside coronal holes, accompanied by low (as compared to 
adjacent quiet sun areas) rate of dipole emergence (Abramenko et al, 
2006; Hagenaar et al, 2008); and ii) magnetic features that contribute 
to the open flux are predominantly well above the noise level, even for 
the noisiest magnetogram. 

For 44 coronal holes analyzed here, magnitudes of the net magnetic 
flux, ^nct = J2B^^As, where As is a pixel size, are presented in the 
5-th column of Table I. The 6-th column shows the densities of the net 
magnetic flux, i?net = J2B\\/A, where A is the area in pixels inside the 
CH boundary, 7-th column in Table I. We would like to emphasize that 
the moduli of Bnet tend to decrease from March 2001 toward July 2006, 
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i.e., toward the end of the 23rd cycle. This is in quahtative agreement 
with Harvey et al. (1982). 

2.4. Correlations Between Calculated Parameters 

Figure 3 shows correlations between parameters of CHs and the solar 
wind speed. The solar wind speed, {Vx), is directly correlated with the 
total net flux (upper left panel). The corresponding Pearson coefHcient 
is 0.71 and the linear fit is given by 

(F,) =494 + 47.4 -^net, (2) 

where (V^) is in km and $net is in 10^^ Mx. 

However, the total net flux is a product of the CH area, A, and the 
average net flux density, i?net- According to Figure 3, the solar wind 
speed is strongly correlated with the CH area (upper right panel, corre- 
lation coefficient (cc)=0.75) and only weekly depends on the averaged 
net flux density (lower left panel, cc=0.20). Linear fitting to the data 
point gave us the following relationship between the speeds and CH 
areas: 

(Vx) = 486 + 8.50 • A, (3) 

where A is in 10^ arcsec^. As a result, there is a very strong dependence 
of the total net magnetic flux from the CH area (cc=0.92, lower right 
panel), which allows a reliable estimation of the total open flux in a 
CH: 

$net = 0.045 + 0.156 • A, (4) 

where $net is in 10^-*^ Mx and A is in 10^ arcsec^. 

These findings very well agree with previously reported results (Wang 
and Sheeley, 1990; Robbins et al, 2006; Vrsnak et al, 2007; Schwadron 
and McComas, 2008) in spite of differences in used techniques. This 
consistency indicates reliability of the coronal hole analysis presented 
here. 

On many occasions, analyzed CHs had closed loop structures embed- 
ded within them, which, in general, reduces the area occupied by the 
open flux and may lead to different magnitudes of the averaged total 
flux density. We excluded the closed loop areas from our calculations 
and repeated the above analysis. As expected, since the closed loop 
regions substitute a small fraction of the CH area, accounting for them 
did not lead to a significant change of the relationship between the solar 
wind speed and the total flux density. In fact, the Pearson coefficient 
even slightly decreased to 0.15. 

We finally would like to note that the well defined relationship be- 
tween the speed and area opens up a possibility to estimate the solar 
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wind speed 3-4 days in advance by measuring the area of a CH, when 
it passes the central meridian. (Note that in this study the CH area 
was measured inside a contour of 80 DN in EIT/Fe XV 284 A images.) 




Figure 3. Correlations between magnetic parameters of coronal holes and solar wind 
parameters. Solid straight lines represent the best linear fit to the data points. 



3. Multifractality in the Magnetic Field in Coronal Holes 

The absence of a relationship between the solar wind speed and the 
density of the magnetic flux, seems to be suspicious at the first sight. 
The solar wind speed represents the intensity of the solar wind ac- 
celeration processes, i.e., energy release dynamics, which ultimately is 
related to the magnetic field. Numerous recent studies of the dynamics 
of energy release in the low corona based on Hinode data (Baker et al, 
2008; Suematsu, 2008; Shimojo, 2008; Moreno-Insertis et al, 2008) tell 
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us that there is a role for the magnetic field to play in these phenomena. 
Two things could be done to resolve the problem. First, improve spatial 
resolution of magnetic field measurements, and second, explore the 
proper measures of complexity of the magnetic field. 

The Solar Optical Telescope (SOT, Tsuneta et al, 2008) onboard 
Hinode has a 50 cm aperture mirror and is the largest optical solar 
telescope ever sent to space. The Hinode/Spectropolarimeter (SP, Ichi- 
moto et al, 2008) obtains the full Stokes parameters using the Fe 
I 630 nm absorption line and offers a unique opportunity to obtain 
magnctograms with pixel size of 0.32 arcsec in the fast mode and 

0. 16 arcsec in the normal mode. Coronal holes at low latitudes are 
rare during a solar minimum, where we are now. Nevertheless, on 11 
November 2008 the SOT/SP instrument recorded a CH (Figure 4) at 
the disk center in the fast mode. We analyzed here one of the SOT/SP 
magnctograms of 942x500 pixels taken in the fast mode (Figure 5). The 
inversion code and calibration routine were performed at HAO CSAC 
{http : / /www. csac.hao.ucar.edu/). The magnetogram was carefully 
despiked. 

It is believed now that magnetic structures in active regions and in 
quiet sun areas arc fractals (Tarbell et al, 1990; Schrijver et al, 1992; 
Balke et al, 1993; Lawrence et al, 1993; Meunier, 1999; Ireland et al, 
2004; Janssen et al., 2003; McAteer et al., 2005). In particular, Meunier 
(1999) and McAteer et al. (2005) showed that magnetic structures of 
active regions, even those recorded with 4 arcsec resolution (MDI full 
disk magnetograms) , are fractals. Janssen et al. (2003) showed that 
the magnetic field in a quiet sun ares is also a fractal. Fractals are 
self-similar, porous structures with jagged boundaries. Their scaling 
parameters do not vary with spatial scale. For example, the filling 
factor, i.e., the ratio of the area occupied by (above-noise) fields to the 
entire area, does not vary when the resolution changes. An example: a 
comparison of the 3rd and 5-th rows in the last column in Table II tell 
us that the fraction of the area occupied by strong fields is nearly the 
same, about 6%, for the 4 arcsec and 1.2 arcsec resolutions. This might 
indicate that the coronal hole magnetic field structure at these scales 
can be considered as a single fractal, or, in other words, a monofractal. 

Temporal variations in monofractal structures are non-intermittent, 

1. e., high fluctuations in energy release are very rare and they do not 
represent a burst-like behavior. Therefore, strong energy release events 
are rare and they do not define the energy balance dynamics. 

Another situation is when a multifractal is formed. Multifractals 
form in nature ubiquitously when several processes contemporaneously 
govern formation of a structure, each one dictating its own rules of 
clustering and fragmentation. A highly intermittent temporal behavior 
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is inherent for multifractals, so that high fluctuations in energy release 
are not rare, and the regime of violent, burst-like energy release dy- 
namics is set. Therefore, to relate the energy release dynamics with 
the magnetic field, one should find the spatial scales where multifractal 
properties of the magnetic field manifest themselves, if any. (Note that 
the terms multifractality and intermittency describe the same property 
of a structure. However, the former is related to the spatial domain 
while the later is usually utilized for the temporal domain; see, e.g., 
Abramenko (2008) for details.) 

One of possible ways to diagnose multifractality is based on calcu- 
lation of the filling factor as a function of spatial scale (Frisch, 1995). 
For multifractals, the filling factor decreases as the scale decreases. In 
other words, the fraction of a volume occupied by strong fields (the so 
called active mode), decreases as we study the multifractal at smaller 
and smaller scales. For monofractals, this ratio is constant with scale. 
To this end, our goal was to explore how the filling factor varies with 
scale inside coronal holes. 

The ratio of the active mode to the entire volume as a function of 
the scale, r, can be calculated from the flatness function, which for the 
longitudinal magnetic fleld Bi, can be written as (Frisch, 1995): 

F{r) = S,{r)/{S2{r)f. (5) 

S,(r) = (|BKx + r)-SKx)n (6) 

are the structure functions, and x is the current pixel on a magne- 

togram, r is the separation vector (i.e., the spatial scale), and q is 
the order of a statistical moment, which takes on real values. The 
angular brackets denote averaging over the magnctogram. Details of 
the calculation routines and applications can be found in Abramenko 
et al. (2002, 2003, 2008) and in Abramenko (2005). The filling factor 
is then calculated as 

/(r) = 1/F(r). (7) 

As it was mentioned above, the filling factor does not depend on the 
spatial scale, r, in case of a monofractal structure. On the contrary, for 
a multifractal, the filling factor displays a power-law decrease as the 
scales become progressively smaller. 

Results of the filling factor are presented in Figure 6. Hinode/SOT/SP 
data for the CH observed on 30 November 2008 (blue curve) indicates 
that at scales larger than approximately 2 Mm, the magnetic field 
structure seems to be a monofractal. Only a very slight slope (0.024) of 
the power law linear fit is observed. Thus, at scales larger than 2 Mm, 
the magnetic field in a coronal hole seems to be a monofractal. 
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Figure 4. The 11:42UT XRT image on 30 November 2008. The box in the cen- 
ter of the solar disk outUnes the area inside a coronal hole where the SOT/SP 
magnetogram was taken (Figure 5). 



For comparison, we calculated and plotted in the same graph (Figure 
6) the filling factor for an active region NOAA 0930 and a plage area 
to the west of this active region. The corresponding magnetograms 
were also derived with the SOT/SP instrument in the fast mode and 
processed with the same routines as those applied for the CH magne- 
togram. Data for the AR and plage show a steeper slope of the power 
law linear fit with indices of 0.18 and 0.09, respectively. The dependence 
of the filling factor on the scale implies multifractality of the magnetic 
field at scales larger than 2 Mm. 

To double check our inference on monofractality in CHs at large 
scales, we calculated the filling factor from the MDI/HR data for 36 
square areas located inside 19 coronal holes observed during 2002-2004 
at the solar disk center. The MDI results are shown in Figure 6 with the 
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Figure 5. Hinode SOT/SP magnetogram of the coronal hole area outlined in Figure 
4 obtained in the fast mode (pixel size of 0.297x0.320 arcsec) and calibrated with the 
HAO inversion code. The magnetogram was despiked. The size of the solar surface 
is 200x115 Mm. The image is scaled in the range from -500 to 500 G. 



dotted line. At scales larger than approximately 3 Mm, the filling factor 
is constant, which confirms our inference on the monofractality of the 
coronal hole magnetic field. At spatial scales smaller than r k, 3Mm the 
filling factor function is influenced by noise and resolution of the MDI 
data and it tends to grow as the scale decreases. This effect is caused 
by the Gaussian nature of the data noise, with its intrinsic filling factor 
of 1/3. 

What is interesting is that the Hinode data do not show similar 
behavior of the filling factor at the same scales, r < 3 Mm. Instead, 
below r = 2 Mm, the Hinode SOT/SP filling factor displays a rapid 
decrease, which does not seem to be caused by data noise. The decrease 
of /(r) is well pronounced for the coronal hole, an active region and 
plage area data. 

The break in the filling factor function is observed at scales ap- 
proximately of 1 Mm. We may conclude that below this scale, the CH 
magnetic field is a multifractal. 



4. Summary and Discussion 

In summary, analyzing magnetic fields inside the low-latitude coronal 
holes, we arrived at the following conclusions. 
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Figure 6. Variation of the filling factor, /, with spatial scale, r. Blue - data from the 
Hinode/SOT/SP magnetogram of the coronal hole shown in Figure 5. Black dots 
- the averaged filling factor calculated from 36 magnetograms for 19 coronal holes 
observed between 2002 and 2004 with SOHO/MDI in the high resolution mode. 
The flat filling factor is observed at scales above 3 Mm in both data sets. The 
decreasing of the filling factor at scales below 2 Mm is observed in SOT/SP data. 
For comparison, data for the active region NOAA 0930 (red) and weak plage area 
(green) derived from the Hinode/SOT/SP fast mode magnetograms are shown. The 
active region and the plage area display the decreasing filling factor at large scales 
above 2 Mm. Dashed lines represent the best linear fit to the data points inside an 
interval of Ar = (1.6 — 8.2) Mm. The steeper slope of the fit corresponds to higher 
degree of multufractality. 

i) The density of the net magnetic flux does not correlate with the 
corresponding in situ solar wind speeds. At the same time, both CH 
area and total net flux correlate very well with the solar wind speed and 
the corresponding spatial Pearson correlation coefficients determined 
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for 44 CHs are 0.75 and 0.71, respectively. The relationship between 
the CH areas and the solar wind speed with almost the same results in 
correlation was earlier derived by Robbins et al. (2006) and by Vrsnak 
et al. (2007). 

As we discussed in Section 1, the fact that the net magnetic flux 
density is not correlated with the solar wind speed can hardly be 
interpreted as irrelevance of the magnetic field to the solar wind ac- 
celeration process. One would rather suggest that the net flux density 
measured with the resolution of 4 arcsec does not reflect the energy 
release dynamics inside CHs. Analysis of SOHO/MDI/high-rcsolution 
and Hinode/SOT/SP magnetograms seems to support this suggestion. 

ii) The filling factor as a measure of multifractality in CHs was cal- 
culated as a function of spatial scale. It was found to be nearly constant 
at scales above 2 Mm. Its magnitude is approximately 0.04 from the 
Hinode data and 0.07 from the MDI/HR data. At scales smaller than 
2 Mm, the filling factor starts to decline as the scale becomes smaller, 
and at approximately 1 Mm the regime of fast decrease of the filling 
factor is set. 

A constant filling factor at r > 2 Mm in the CH magnetic fields indi- 
cates their monofractal nature and self-similarity. A self-similar struc- 
ture by definition possesses constant statistical parameters at all scales 
(such as various scaling exponents, including the filling factor). How- 
ever, only for artificial, mathematically created fractals, self-similarity 
is present at infinite range of scales (see, e.g., Schroeder, 1991). For 
a majority of structures in nature, however, self-similarity with con- 
stant scaling parameters only holds at a finite range of scales, while 
at the entire interval of scales the scaling parameters are different. As 
a result, a multifractal structure forms with a superposition of many 
fractals, each one imposing its own scaling rules. A crucial difference 
between monofractals and multifractals is in their temporal evolution. 
In monofractals, large fluctuations of parameters (say, energy release 
events) are rare and do not determine mean values. In other words, time 
profiles are non-intermittent and evolution proceeds without catastro- 
phes. On the contrary, in multifractals, the time profiles are highly 
intermittent, large fluctuations are not rare, and they determine mean 
values. The temporal energy release process is burst-like. 

If so, the monofractal property of the CH magnetic field at scales 
above 2 Mm seems to be the most plausible reason why the averaged 
magnetic flux density, derived from observations with low resolution, 
does not correlate with the solar wind speed: the bulk of energy release 
dynamics, needed for the solar wind acceleration, occurs at smaller 
scales, where the magnetic field structure is entirely different. 
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Berger and colleagues (Berger et al, 2004) observed a plage area 
with 0.1 arcsec resolution, and they report that the magnetic flux is 
structured into amorphous ribbon-like clusters with embedded knots 
of enhanced density, which seems to correspond to the notion of mul- 
tifractality at small scales. 

It is worth to mention that the property of monofractality of solar 
magnetic fields was known for a long time (see references in the previous 
section). Difficulties of describing magnetic structures with a single 
fractal dimension at the entire range of scales was also noticed (see, e.g., 
Tarbell et al, 1990; Janssen et al, 2003), which actually is a signature of 
multifractality. Fractal analysis of high resolution magnetograms from 
VTT with 0.4 arcsec spatial resolution for a quiet sun area (Janssen et 
al, 2003) revealed a break of self-similarity at scales of 1.3 Mm, which 
is very close to the scale found in this study. 

For an active region and a plage area, our approach for deriving the 
filling factor of the magnetic field produced that /(r) = 0.14 — 0.17 
at r = 0.3 Mm that generally agree with earlier reports. According to 
Tarbeh et al. (1979), Schrijver (1987), Berger et al. (1995) and refer- 
ences herein, the magnetic filling factor is typically inside a range of 10 
- 25%. In active regions, where the range of solar flares spreads over all 
observable scales, multifractality is also present at the same range of 
scales. This presents further evidence that energy release dynamics and 
the multifractality are mutually related properties of solar magnetic 
fields. 

With the Hinode instrument in operation, many new phenomena 
related to the CH dynamics, coronal heating and solar wind acceler- 
ation will be discovered now, when the spatial scales less than 1000 
km are available for analysis. Recent study of the evolution of network 
magnetic elements (Lamb et al, 2008) proved that processes at sub- 
resolution scales are of vital importance for understanding the observed 
dynamics of magnetic flux. 
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CH 
name 


Culmination 
time 


^0 5 

dov 


tA, 

doy 


$net, 

10^^ Mx 


B 


G 


A, 
pixels 


deg 


km 


CH191 


2001Mar02/23:00 


61.96 


64, 


.96 


3. 


,00 


5. 


,01 


29145 


17, 


,2 


559 


±16 


CH230 


2002Feb03/22:00 


34.92 


37 


.75 


1. 


,09 


3. 


.43 


15408 


-2. 


,8 


614 


±27 


CH201 


2002Mar02/17:00 


61.71 


64, 


.38 


3. 


.03 


4. 


.27 


34472 


2. 


.2 


649 


±28 


CH202 


2002Mar29/ 18:00 


88.75 


91, 


.08 


2. 


.61 


4. 


.22 


30127 


1. 


.7 


701 


±38 


CH210 


2002Apr28/23:00 


118.96 


122, 


.63 


0. 


,84 


3, 


,85 


10621 


-19, 


,6 


462 


±14 


CH231 


2002Jul03/12:00 


184.5 


187, 


.67 


-0, 


.92 


-4, 


.60 


9721 


-13. 


,2 


536 


±16 


CH204 


2002Nov01/17:00 


305.71 


309 


.04 


1. 


.49 


2. 


.54 


28578 


-6. 


.3 


549 


±16 


CH350 


2003Jan27/06:00 


27.25 


30, 


.67 


-0, 


.24 


-4, 


.75 


2555 


1, 


,6 


507 


±14 


CH302 


2003Feb23/19:00 


54.79 


58, 


.13 


-2. 


.04 


-4. 


.26 


23176 


7. 


.1 


560 


±18 


CH303 


2003Mar01/03:00 


60.13 


63, 


.29 


-1. 


.65 


-4. 


.76 


16856 


14. 


.2 


530 


±32 


CH304 


2003Apr08/09:00 


98.38 


100, 


.79 


3. 


.12 


2. 


.73 


55633 


-7. 


.9 


702 


±27 


CH305 


2003Apr24/07:00 


114.29 


117, 


.71 


-0, 


.29 


-2, 


.67 


5407 


0. 


,2 


510 


±23 


CH307 


2003May03/19:00 


123.79 


126, 


.38 


3. 


,07 


3. 


,23 


46245 


-11. 


,1 


693 


±23 


CH354 


2003May21/01:00 


141.04 


144, 


.54 


-1. 


.43 


-3. 


.12 


22381 


-1. 


.0 


530 


±31 


CH309 


2003May25/19:00 


145.79 


148, 


.21 


-0. 


.93 


-5. 


.59 


8151 


-20 


701 


±31 


CH312 


2003Junll/ll:00 


162.46 


165, 


.63 


-0, 


.47 


-2, 


.51 


9260 


-0, 


,6 


548 


±21 


CH355 


2003AuglO/01:00 


222.04 


224, 


.63 


-2, 


.61 


-2, 


.43 


52146 


3. 


,6 


650 


±25 


CH314 


2003Aug20/19:00 


232.79 


235 


.04 


4. 


.57 


2. 


.66 


83592 


-4. 


.9 


741 


±28 


CH356 


2003Aug30/03:00 


242.13 


245, 


.29 


-0, 


.40 


-3, 


.75 


5191 


7, 


,8 


526 


±20 


CH357 


2003Sep07/23:00 


250.96 


253, 


.63 


-3, 


.17 


-4, 


.60 


33515 


-19. 


.2 


638 


±28 


CH315 


2003Sepl6/08:00 


259.33 


261.5 


4. 


.84 


2. 


.98 


79079 


-4. 


.2 


749 


±35 


CH317 


2003Dec01/05:00 


335.21 


339, 


.21 


0. 


.12 


1. 


.90 


3207 


-6. 


.0 


449 


±22 


CH320 


2003Decl8/19:00 


352.79 


355, 


.79 


-4, 


.49 


-3, 


.98 


54777 


6. 


,3 


610 


±19 


CH455 


2004Mayl7/10:00 


138.42 


141, 


.75 


0. 


,21 


2. 


,07 


5074 


-2. 


,7 


519 


±24 


CH456 


2004May28/23:00 


149.96 


153, 


.13 


-0. 


.58 


-1. 


.17 


24452 


13. 


.0 


549 


±17 


CH457 


2004May31/03:00 


152.13 


155, 


.54 


-3. 


.03 


-4. 


.24 


34826 


-3. 


.3 


503 


±15 


CH458 


2004Jun05/01:00 


157.04 


160 


.96 


-0, 


.40 


-3, 


.73 


5281 


14, 


,2 


452 


±13 


CH459 


2004Junl2/10:00 


164.42 


167, 


.75 


0. 


,31 


1. 


,83 


8469 


-9. 


,8 


542 


±25 


CH460 


2004Jul07/23:00 


189.96 


193 


.71 


0. 


.36 


2. 


.36 


7540 


11. 


.4 


491 


±20 


CH461 


2004Jull3/13:00 


195.54 


199 


.13 


0, 


,15 


1, 


,32 


5590 


2, 


,8 


518 


±22 


CH462 


2004Nov21/19:00 


326.79 


330 


.21 


0. 


.28 


3. 


.53 


3967 


16. 


.0 


513 


±17 


CH412 


2004Nov27/17:00 


332.71 


335, 


.63 


-2. 


.47 


-4. 


.29 


27994 


-1. 


.2 


638 


±19 


CH502 


2005Febl5/07:00 


46.29 


49, 


.54 


-0. 


.80 


-3. 


.24 


12051 


14. 


.8 


547 


±21 


CH503 


2005Aprl7/02:00 


107.08 


110, 


.42 


0. 


,45 


2, 


,90 


7650 


-9. 


,5 


523 


±22 




zuuoiviayzy/ uo.uu 




153 


.63 


0. 


,50 


z. 


Oft 




-14. 


,0 


434 


±16 


CH558 


2005Jul26/10:00 


207.42 


210, 


.58 


-0. 


.78 


-2. 


.08 


18438 


-5. 


.4 


559 


±31 


CH563 


2005Octll/19:00 


284.79 


289, 


.63 


0. 


.13 


1. 


.30 


5075 


1. 


.9 


372 


±19 


CH564 


2005Oct23/15:00 


296.63 


299 


.88 


0. 


,24 


1, 


,83 


6557 


-1. 


,9 


508 


±26 


CH567 


2005Dec26/15:00 


360.63 


363 


.04 


-0, 


.70 


-1, 


.75 


19465 


2. 


,3 


668 


±32 


CH601 


2006Janl2/19:00 


12.79 


16 


.71 


0. 


.41 


1. 


.95 


10396 


4. 


.3 


433 


±12 


CH603 


2006Febl2/03:00 


43.13 


46 


.79 


0, 


,29 


2. 


,48 


5799 


-8, 


,3 


535 


±21 


CH604 


2006Feb24/19:00 


55.79 


60 


.21 


0. 


.40 


4. 


.14 


4813 


-4. 


.9 


411 


±14 


CH606 


2006May03/09:00 


123.38 


126, 


.88 


0. 


.78 


1. 


.33 


28904 


-4. 


.0 


600 


±27 


CH607 


2006Jul02/19:00 


183.83 


186, 


.67 


-0. 


.91 


-1. 


.30 


34120 


2. 


.0 


593 


±24 
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Table II. SOHO/MDI magnetic measurements for a box inside CH204 coronal hole. 



Type of 


-Bnet 


Btot 




Flux 


Bnet{<p) 


Area above 


magnetogram 


G 


G 


G 


imbalance, % 


G 


noise level, % 


Full disk 1 min 


4.53 


19.5 


17.3 


23 


1.11 


4.1 


Full disk 1 min 














smoothed 


4.50 


10.3 


4.17 


44 


0.38 


19 


Full disk 5 min ave 


4.48 


12.1 


9.3 


37 


0.84 


6.1 


Full disk 5 min ave 














smoothed 


4.48 


8.1 


3.2 


55 


0.23 


16 


Hi res 


3.88 


12.4 


8.9 


31 


1.64 


6.6 
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